Induced superconductivity in 2D electronic systems 
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The approach applicable for spatially inhomogeneous and time-dependent problems associated 
with the induced superconductivity in low dimensional electronic systems is developed. This ap- 
proach is based on the Fano- Anderson model which describes the decay of a resonance state coupled 
to a continuum. We consider two types of junctions made of a ballistic 2D electron gas placed in 
a tunnel finite-length contact with a bulk superconducting leads. We calculate the spectrum of the 
bound states, supercurrent, and the current-voltage curve which show a rich structure due to the 
presence of induced gap and dimensional quantization. 
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I. INTRODUCTION 

Recent progress in studies of transport in graphene 
that followed the seminal work of Refi has boosted the 
interest in properties of contacts between the graphene 
sheets and various types of electrodes attached to them. 
Of special interest are contacts between graphene and 
superconducting electrodes due to the specific nature of 
the Andreev reflection in graphene^. In practical devices, 
such type of contacts has the form of a superconduct- 
ing lead placed on top of the graphene layer which par- 
tially overlaps with the lead. Following Refi^, the model 
most commonly used for description of such contacts as- 
sumes that the superconductor simply introduces certain 
pairing potential in the part of layer which is immedi- 
ately under the superconductor as well as shifts its Fermi 
level away from the Dirac point. As a result, the con- 
tact is treated as being formed between the usual normal 
graphene layer and such piece of graphene where both 
the induced pairing potential and the high doping level 
are present. Various modifications of this model have 
been extensively used for studying transport properties 
of graphene contacts 3 - - — . 

Though this model is a significant step forward in un- 
derstanding the induced superconductivity in graphene, 
it still oversimplifies the proximity effects which the bulk 
superconductor has on the underlying normal sheet. To 
study these effects more carefully one can look at other 
models used for contacts between superconductors and 
low-dimensional electronic systems. For example, in 
Refsj 10 ! 11 the proximity effect in a two-dimensional (2D) 
electron gas was modelled by a uniform plane contact 
between a bulk superconductor and a thin normal con- 
ducting layer. Such model correctly catches the main 
physics of the proximity effect in spatially homogeneous 
structures in equilibrium. In particular, it accounts for 
the induced energy gap that depends on the contact resis- 
tance. However, it cannot describe spatially-dependent 
problems in 2D systems nor can it be easily extended to 
time-dependent or non-equilibrium phenomena. 

In the present paper we develop an approach which 



is suitable for spatially inhomogeneous and/or time- 
dependent problems associated with the induced super- 
conductivity in 2D systems placed in a contact with a 
bulk superconductor. Our model is similar to that used in 
RefSi 12 : 13 for impurities in a superconductor and is based 
on the so-called Fano-Anderson model which describes 
the decay of a localized state coupled to a continuum 14 . 
Our model can be applied to various 2D electronic sys- 
tems, including simple 2D gas, graphene layer, etc. In 
this paper we consider two particular examples of junc- 
tions made of a 2D ballistic electron gas placed under 
the superconducting electrodes. The application of our 
model to the induced superconductivity in graphene will 
be considered elsewhere. 

The paper is organized as follows. In the next Section 
we describe our proximity model in its general formu- 
lation suitable for various applications. The particular 
example of a junction between two proximity-induced 
superconductors is considered in Section IIII1 We solve 
equations for the Green functions, find the energies of 
the bound states, and calculate the supercurrent through 
such junction. In Section IIVI we calculate the current- 
voltage curve for a contact between the semi-infinite nor- 
mal and the proximity-induced superconducting regions. 
Our results are summarized in Section fVl 



II. MODEL 

Consider a 2D electron layer placed under a bulk su- 
perconductor and coupled to it via tunnelling through a 
thin insulator coating, Fig. [T] The Hamiltonian of the 
system has the form 



H — Hs + H2D + Ht 
In the superconductor 



(1) 



He 
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J2*U*Ws-E F ]i> a (r) 



AM(r)ft(r) + A*^(r)* t (r) 
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FIG. 1: Superconducting correlations in the 2D system under 
the superconductor are induced through a thin insulator. 



where Ef is the chemical potential in the superconduc- 
tor, e s is the kinetic energy operator. For parabolic spec- 



The coordinate r is 



trum it is e s = ^ 
a 3D vector which belongs to the superconductor. In the 
2D layer 



H 2D = dJ d 2 R Si (R) [c2d(R) - E F ] a a (R) , (3) 



where d is the layer thickness, e 2 £)(R) is the kinetic en- 
ergy operator in the 2D layer. The coordinate R is a 2D 
vector which belongs to the layer. The creation and anni- 
hilation operators in the layer d^ a (R) , a a (R) are normal- 
ized to the layer volume, such that the anti-commutator 
[a a (K),al (R')] + = d _1 5(R-R'). The tunnel Hamilto- 
nian has the form 



H T =dY / J [*l,(r)T(r,R)a a (R) 



it(R)Tt(R,r)* Q (r) 



d 3 rd 2 R 



(4) 



Here the coordinate R refers to the layer while coordi- 
nate r refers to the superconductor; the matrix element 
T(r, R) describes tunnelling between the layer and the 
superconductor, T*(r, R) = T^(R, r). This model is sim- 
ilar to the model used in Refs i 12 ' 13 for an impurity in 
superconductor . 

The Matsubara Green functions are 



T T d a (R 1 )dt(R 2 ) 



<5q^G2_d(Ri, R2) , 



T T tf Q (n)a£(R2)) =<5 a/3 G T (ri,R 2 ) , 



T T # Q (n)M(r 2 )) = 5 a pG s (rx,r 2 ) , 



and 



(T T a Q (R 1 )a (3 (R 2 )) = ia i ^F 2D (R 1 ,R 2 ) 



r T * Q (n)o^(R2) 



T T * Q (n)^(r 2 )\ =ia^F s (ri,r 2 ) , 



etc. We introduce the Nambu matrixes 



Hs(r) 



e s - E F -A 
A* e s - E F 



, G = 



G F 
-Ft G 



and denote 

GsHr)=f 3 h-^- + Hs(v) 

the inverse operator in the superconductor. If needed, it 
can also include the impurity self energy. Neglecting the 
back-action of the thin 2D layer on the superconductor, 
we have for the superconducting Green function 

<% x (ri)(5 s (ri,r2) = lM(n - r 2 )rj( Tl - r 2 ) . 

Equations for the mixed Green functions Gt (i"i,R 2 ) 
can be written in the form 



G^ 1 (r 1 )G T (r 1 ,R 2 ) + dy f(r u R')G 2D (R', R 2 ) d 2 R' = 
where 

T(r,R) 



T(r,R) 







T*(r,R) 
This gives 
G T (ri,R 2 ) 

= -df G s (r 1 ,r')f(r',R')G 2D (R',R 2 )d 2 R'd 3 r' . (5) 



Equations for the Green functions in the layer can be 
written as 

G 2 ^(Ri)G 2D (Ri,R 2 ) + / f\R 1 ,r')G T (r',R 2 )d 3 r' 
= ld-^^Ri - R 2 )5(n - r 2 ) 



where the inverse operator in the layer is introduced, 

G-^R) =m^ + MB) - E F ] . 
Here the kinetic energy operator is 



1 ( - t 9 - e * 1 
e 2 D = ^— I -in,—— - r 3 -A + e , 



(6) 



2m 



9R 



while eo is the shift of the bottom of the 2D conduction 
band measured from that in the superconductor. 
Using Eq. ([5]) this equation becomes 

G 2 d(Ri)G ! 2 d(Ri,R 2 ) — J St(Ri,R') 
xG 2 d(R', R 2 ) d 2 R' = 1M- 1 ,5(R 1 - R 2 )6(t! - t 2 ) (7) 
where 

t T (R u R')=df f^(R 1 ,r')G s (r',r")f(r",R')d 3 r'd 3 r". 

We assume a site-to-site tunnelling which does not con- 
serve momentum, T(r,R) = £<5(r 11 — R)<5(z), etc., where 
i is real and does not depend on R. Here 8(z) selects an 
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average value of a function at a distance of the order of 
inter-atomic scale near the surface. As a result, 

yi /ft r> G s (R,R,;0) F s (R,R,;0) 

L T (R i) R J -)-*^_ j pt (R . !R .. 0) ^(R^R^O) 

(8) 

where (Rjj 0) means that r\\ = Ri and — z,j = 0. 

The coordinate dependence of the Green functions has 
the formic 



of the Matsubara functions onto the real-frequency axis. 
We introduce the Keldysh matrixes 



fin r<K 

Ho %* 1 5 



t A 



G s (r t ,r j ) = 



2irh?\ri — Tj\ 



[lcos(pir|rj - rj\/h) 



The equations become (the index 2D is dropped) 

(e 2D -E F - f 3 e)e e>e ,(R,R') - [<S(R)S(R,R')] e>e , 

= itfd^SCR - R)2irS(e - e') (12) 

where 



+ ism( PF \r t -r J \/h)g s (h,r)} . (9) [S(R)0(R, R% e , = [ S e>ei (R)& li6 ,(R, R')|% 



Here £ is the mean free path is the superconductor, n = 
(r.; — rj)/\(ri — Tj)\ is a unit vector, r = (r^ + fj)/2, and 
gs is the quasiclassical superconducting Green function 
integrated over the energy variable, Eq. (fTOl) . 



<7s(PF,r) = / ^G s (p,r) 



(10) 



which depends only on the direction of the momentum 
p F = p F n. 

The self-energy Ex(Ri,R;) in Eq. (|SJ oscillates as 
exp(ipi?|Rj — Ri\/h) while the Green function contains 
exp[ip(R;-R^.)/ft] = exp[ip(Ri-Rj)/h+ip(Rj- R'^/h] 
where |p| is close to the Fermi momentum p 2 u m the 2D 
layer. If p2D ^ p F , the integral in Eq. ([7]) converges at 



R — Ri — R; 



where a is the interatomic distance. 



Therefore, one can invoke the approximation 

E T (R„R) = S(R,)<5(R J -R) . 

The first term in Eq. (j9]) diverges for |R— R'| — > 0, but it 
does not depend on the superconducting state. It is real 
and has the form of an effective potential. We can include 
it into the shift of the bottom of the conduction band eo 
or into the chemical potential assuming that /i already 
accounts for it. The remaining term for |R — R'| — » 
gives 



E(R) = »r (5s(R;0)) 



(ii) 



Here we introduce the tunnelling rate T = irv^dsot 2 
where = mp F /2TT 2 h 3 is the 3D density of states, 
so ~ o 2 is of the order of the area of the unit cell in the 
layer; the overall order of magnitude is V ~ t 2 /E F . An- 
gular brackets denote averaging over momentum direc- 
tions. The quasiclassical Green function gs(R;0) taken 
at the superconductor/2D-layer interface is the only char- 
acteristic which is needed in our model to account of the 
properties of the bulk superconductor. Since the qua- 
siclassical function varies over distances of the order of 
superconducting coherence length, the exact atomic-scale 
boundary conditions at the interface for the microscopic 
wave functions are not critical for our model. 

We write the final equation in 2D layer for the real- 
time Green functions making the analytical continuation 



III. S/2D/S STRUCTURE 

A. Green functions 

In this paper we restrict ourselves to stationary prob- 
lems and put G f[f^ = G^ A ^2irH5(e— e'). In the following 
Sections we apply our model to the S/2D/S structures 
which are made of a ballistic 2D electron gas placed in a 
contact with bulk superconducting leads (S). One of the 
examples is shown in Fig. [2] A generic structure of this 
type consists of a layer of 2D electron gas which is infi- 
nite in y direction and has boundaries at certain x where 
R = (x, y). Transforming to the plane waves along the y 
direction, 

Gf(RR') = J ^(y-v')/K G K py {x,x>)^ , 

we observe that, as functions of x for x ^ x 1 , the re- 
tarded Green functions G R Py (x,x') and (x, x') sat- 
isfy the set of 2 linear homogeneous second-order differ- 
ential equations which follow from Eq. (fT2l) . 

t^d^_ 

2m dx 2 

2m dx 2 



u(x) - (e + rji)u(x) + ri 2 v(x) = 0, (13) 
v(x) - (e + rji)v{x) + n\u{x) = 0, (14) 



where fi = E F — eo is the chemical potential in the 2D 
gas measured from the bottom of the conduction band, 
Mx = I 1 ~ Py/2m, and 

m = (g R ) , m = ^^ </«) ,4= t T {ft R ) . (15) 



s 


2D 


S 









-(L+D/2) 



FIG. 2: Contact made of two superconducting electrodes 
placed on top of the 2D electron gas at a distance D from 
each other. Each electrode covers a length L of the 2D layer. 
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Eqs. (|T3|) . (|14l) have 4 independent solutions, 

i = 1,2,3,4 



Using the standard method in the theory of differential 
equations, one can write the Green functions as sums 

( F*(x?J) ) = CiW&VW + C 2 (x')& 2 Hx), x > x', 
( F^x*')) ) = ~ C 4 (x')* (4) (x), x < x'. 

The functions are chosen such that the retarded 

Green functions are regular in the upper half-plane, 
Ime > 0. The coefficients Ck(x') do not depend on x 
but can be functions of x'. We find from Eq. (IT2l 

G(x + 0,x)-G(x-0,x) = , 

* G ( X + 0iX) -JL G (x -0,1) = -^, 
ax ax na 

F^(x + 0,x) - F*(x-0,x) = , 

4^ t (a; + 0,x)--^-F t (a:-0, : E)=0. 
ax ax 

Let us define the vectors 



Ci = 



and the matrix \J\ where 




Rk = — TT 




C/j = U W , = dt»W /da: , f/j = , = dt> (l) /da . 

The upper index numerates columns, while the lower nu- 
merates rows. Equations take the form J^i^ld = R-k- 
The solution is 



k 



(16) 



Here W — det [Z7£] is the Wronskian and A\ (x) = 
{-l) i+k W l k {x) where W l k {x) is a minor of det [£/£]. The 
Wronskian is independent of coordinates. 



Here k± = k x ± e/ , while fc^ and v x are the a; com- 
ponents of the Fermi wave vector and Fermi velocity 
hk x — mvx in the 2D gas, h 2 k 2 /2m — [j, x ; the particle 
and hole Nambu vectors are 



The set of two other functions is obtained using reflection 
from the left boundary, 



L -ik-(x+B-) 



(20) 



In the right induced-superconductivity region D/2 < x < 
L + D/2, the four functions satisfying Eqs. (JT3J) - (fi"4|) 
are linear combinations of 




(x — L — ^- ) 




The functions in the left superconducting region are ob- 
tained by replacing <fi with —<j> and {L + D/2) with 
— (L + D/2). Note that the momentum in the induced- 
superconductivity region, 



k s ± = k x ± j^-\/'< - '/i y 



(21) 



has the same Fermi-momentum projection on the x coor- 
dinate as in the normal region. Thus the normal reflec- 
tion at x = ±D/2 does not occur. The coherence factors 
are 



u(e) - 4=J 1 + — T— 
V2\ e + Vi 




(22) 



B. Reflection coefficients 

Consider structure shown in Fig. [5J The right su- 
perconductor has the phase <j>/2 while the left one has 
—(f>/2. The four independent functions needed for the 
Green functions can be constructed as follows. The first 
two wave functions inside the 2D layer —D/2 < x < D/2 
contain waves incident from the left and reflected from 
the right boundary, 



+f R e ik_ {x -S.)^ 



h ■ 



(18) 



Here 



e = (e + m f - mv \ = 6 2 -T 2 



since rj 2 — 172% = ~T 2 . Moreover 



Ve 2 - A 2 



A 



I 2 = e 2 - 



(23) 



2;r e 2 



Ve 2 - |A| 2 



The square roots are defined as analytic functions in the 
plane of complex e with cuts from |A| to 00 and from 
—00 to — |A| taken at the upper bank of the cut for e > 
|A|. The usual BCS formulas are recovered if one uses 
e — > e + r/i , A — > 779 , A* — > ry 2 such that the BCS gap 
A| 2 is replaced with r 2 |A| 2 /(|A| 2 - e 2 ). 



5 



The boundary conditions require that i*W — = 
at x — ±(L + D/2) together with continuity of 
uW )t ,W , du^/dx , dv^/dx at x = ±D/2. This yields 



i N 

'■a 



N 



r N e lS , f* = rjve" 15 , 



r A e 



■4/2 



' 1 A 



r A e 



i<p/2 



(24) 
(25) 
(26) 
(27) 



where 



r N -- 



( u 2 _ ^ e i(6 + -S-) 
' u 2 _ v 2 e 2i(S+-S-) 



r A - 



u v 



[1 - e»( a +- 4 -)] 



Here we denote <5± = fcJL and (5 = 5 + + 6- = 2k x L, 
while 



5+ - 6- = {2L/hv x )^i 2 -Y 2 



(28) 



One can check that the Green functions defined according 
to Eq. (QH) with the basis functions Eqs. P7 |) -([20 | are 
regular in the upper half-plane of complex e. For e < A 
the coefficients satisfy 



M 2 + \r A \ 2 = 1 , r* N r A + r* A r N = 



(29) 



because there is no quasiparticle flux into the supercon- 
ductor. 

As in RefJ^ one can define an equivalent barrier hight 
Z associated with the end of the conduction layer, 



z/x/i + z 2 



A8+S-) 



(30) 



such that the coefficients have their standard for m 17 ' 18 



r A -- 



u 2 + (u 2 — v 2 )Z 2 



(u 2 -v 2 )zVTTz^ 

TN = 2 i ^ 2 2W2~' ( 31 ) 

u z + {u z — v z )Z z 



The limit L = corresponds to strong normal reflec- 
tion, Z — > oo and |rjy| = 1, while ta = 0. With 
this definition of effective barriers, the problems of the 
bound-state spectrum and of the supercurrent in many 
respects reduce to the corresponding problems in double- 
barrier superconductor/normal/superconductor (SINIS) 
structures. 

Let us denote e g < A the energy at which the square 
root in Eq. (|2"TT) vanishes, i.e., at which £ defined by Eq. 
([23)1 turns to zero, 



(l + 2iy 



A 2 



T 2 = . 



(32) 



The energy e g is the induced gap in the 2D layer. For Y <C 
A we have e g = Y; if Y > A one has e g = A(1-2A 2 /T 2 ). 
The induced gap as a function of Y is shown in Fig. [3] 
For a 2D gas in a direct contact with a superconductor, 
the induced gap was obtained in Ref^. 

If e 2 < e 2 , both the square root in Eq. (|2"Tj) and £ are 
imaginary. The value Z is real and satisfies 

Z/y/l + Z 2 = cxp \-2(L/hv x ) ^/Y 2 - l 2 \ . (33) 



For large L we have Z — 0; the particles do not feel the 
dead end of the superconducting region, such that the 
normal reflection is absent, while \r A \ = 1. This limit is 
realized when L 3> £20 where 

&d = hv 2 D/e g 

is the coherence length of the induced superconductivity, 
V2D is the Fermi velocity in the 2D system v 2 D /2m = [i. 
The limit of long leads is rather hard to realize because 
^2D is considerably longer than £0 in a bulk superconduc- 
tor. On the contrary, for short leads, L — 0, there exists 
only normal reflection, Z — > 00, and |r/v| = 1. 

If e g < e < |A| both the phases 6± and the coherence 
factors u, v are real, while Z is complex. The Andreev 
reflection vanishes while |r/v| = 1 if <5 + — 5— = nn. 

For e > |A| the coefficients still have the form of Eq. 
(f3"Tj) with a complex effective barrier strength Z. For 
large e 3> |A| one has e = e + iY. The effective barrier 
height saturates at Z 7^ 0, while v — > and thus ta — > 0. 
At the same time 



exp [— ALY/hv x ] 



Equation (f2T)|) no longer holds because the quasiparticles 
escape into the superconductor. 



C. Bound states. Supercurrent 

Using the functions Eqs. (JI7J - $U\) and Eqs. ([2"3]l- 
(|2~T1) . we find the Wronskian 

W = 16fc+fc_e 2t7 [sin 2 (/3 + 7 )- |r A | 2 cos 2 (</>/2) 
-\r N \ 2 sm 2 (a + 5)] 

where a = k x D, j3 — eD/hv x , while 7 is the phase of 
the reflection coefficient r^, i.e., e 2 * 7 = rpj/r* N . The 
spectrum is determined by W = 0, which gives 

sin 2 ( / 3 + 7 ) = |r A | 2 cos 2 (0/2) + |rAr| 2 sin 2 a' (34) 

where a' — a + 5. This agrees with the previous cal- 
culations of Ref. 18 for ballistic SINIS contacts. Here we 
concentrate on short contacts such that /3 <C 1. Neglect- 
ing j3 and keeping L finite implies that one has to assume 

If the energy is above the induced gap, e g < e < A, 
the combination 5+ — <5_, as well as u and v are real. 




FIG. 3: Induced gap as a function of the tunneling rate. 
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For energies below the induced gap, e < e g , the barrier 
strength Z, Eq. (I3TJ1) . is real, while u = v* . Therefore, 
Eq. (f34|) yields for the bound state energy 



(e + Vi) 2 _ 1 sin 2 (0/2) 



V2V2 



1 + A 



(35) 



where 



A = 
A = 



sin 2 a' 



; = , e„ < e < A (36) 

sm 2 [2Ly/P-T 2 /hv x ] 



sin 2 a' 



sinh 2 [2Wr 2 - e 2 /hv x ] 



(37) 



Equation (j3"6")l goes into Eq. (j3"?| if one continues ana- 
lytically the phase difference S + — <5_ as a function of e 
around the induced gap e g . 

For energies below A qusiparticles cannot escape into 
the superconductors and form the bound states as a 
result of two processes: the Andreev reflection at the 
boundary between the normal region and the region with 
induced-superconductivity plus the normal reflection at 
the ends of the 2D layer. The spectrum is determined by 
Eq. (|35[): it consists of a series of levels for each p y and 
4> spaced with a distance ~ hv x /L and filling the inter- 
val < e < A. The lowest energy level lies below the 
induced gap if (2LT/hv x ) sin(0/2) > sin a'. For short 
leads, 2LT/hv x -C sin a', the energies are close to the 
levels of geometrical quantization in a potential well of 
length d + 2L, i.e., they satisfy 5+ — <5_ ± a' = xn. If 
L <C hv x /A, there is only one bound state with energy 
e < A. 

The supercurrent is^S 



I = / cos 9 d9 -rv tanh — ^- . 



(38) 



Here 



R 1 = 2e 2 wv 2 V2D/^ = G a wp2D/Tth 



is the Sharvin conductance of the contact, Go = e 2 /irh 
is the conductance quantum, p 2 rj — mV2D is the Fermi 
momentum in the layer, and v 2 = m/2irh 2 is the normal 
density of states in the 2D gas. The conductance is pro- 
portional to the number of modes in the 2D channel of 
width w in the y direction. The sum collects all bound 
states with energies e < A for given p y = p 2 D sin 9. The 
continuum states with e > A give the contribution pro- 
portional to the length D and can thus be neglected. 



1. Short leads 

Consider the case of short leads, LT/hv x -C 1 for 
r <C A. This is the most practical situation in the exper- 
imental devices. If Hv x / A «L« hv x /T, there is a series 
of bound states with energies < e n < A. As discussed 



above, they are close to the levels of dimens ional quan- 
tization. Putting x = (2LT/hv x )y / e 2 /T 2 - 1 we observe 
that x — x n + 5x n where 



sin x n = sin 2 a' , 



and 



Sx n = 



V hv x J 



sin^/2)- 



2x 2 



(39) 



(40) 



which follows from Eqs. (|3"5)) , (|3"6"1) . Equation gives 



n= 1,2,.. 
n = . 



(41) 



Here xq = arcsin | sin a' \ is the smallest root, < xq < 
7r/2, of Eq. (|3T))) . This solution holds as long as Sx n -C 
x n , i.e., for all states with n ^ 0. It also holds for the 
lowest energy state, as long as the level is not very close 
to the resonance, |sina'| 3> (2LT/hv x ). 

Close to the resonance, |sina'| -C 1, the lowest level 
can also be obtained by expanding the sine functions in 
Eqs. ([36 ]l . (|37 |l . Putting i m e we find 



e 2 = 



A 2 + r 2 (i + a 2 



[A 2 -T 2 (l + a 2 



+ A 2 r 2 sin 2 ((/>/2) , (42) 



where a 2 = h 2 v 2 sin 2 a' / (2 LT) 2 . Close to the resonance, 
aT < A, 



Ty/a 2 + cos 2 (0/2) 



(43) 



This goes over into Eqs. (|1D|) when a > 1. The 

typical spectrum is shown in Fig. Ufa). 

If the leads are very short, L <C hv x /A then there 
exists only one level. Except for a very narrow vicinity 
of resonance, one has a 2 T 2 3> A 2 , thus the level is close 
to the superconducting gap A, 



A 2 



(A 2 /a 2 )sin 2 (0/2) . 



(44) 



To calculate the supercurrent Eq. (|38|) we integrate 
the contribution from each level over the incident angle 
9 defined according to v x — v 2 dcos9. Since the reso- 
nance form of the level, Eq. (1431) . holds only in a very 
narrow region of angles, its contribution to the current 
is small. The current is thus mostly determined by the 
off-resonance levels, Eqs. (|3Tj). (|40|. or (J44J) with e„ > T. 
For hv x / A «L« hv x /T, there is a series of levels de- 
fined by Eqs. {H]). ([40]). The current is 



Lr 2 



tt/2 



tan xo 



eRo \hv 2 D 

00 

x F(x ) + ^[F(7rn + a;o) -F(im-x )] 



d9 (45) 
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where 

. i . hv x x 

i^) = _tanh— . 

Due to a large argument in sin(cv') = sm.[k2D(d + 
2L) cos 0], the level xo oscillates rapidly assuming val- 
ues within the interval < xo < ir/2 many times as 
the incident angle varies from to n/2. Let us con- 
sider a function f(x o ,0) of the rapidly oscillating vari- 
able Xq(0) and a slow variable and define the aver- 
age function (/(#)) 80 = J„ + 6 f(xo,0) dO where the in- 
tegral is taken over the full variation range of xq putting 
dxo = (d + 2L)k2D s'm0d0. Since the range < x < n/2 
corresponds to a small variation 60 <C 1, one can keep 
the slow variable constant during integration over dxo , 

(W)) = - I"'* f(xo,0)dx a . 



The integral of the rapid function f(xo,0) can now 
be replaced with the integral of the average function 



For low temperatures, T <C Hv2d/L one can simplify 
the expression f)45f) for the current. We note that for 
n / the integral in (tanxo^X 7 ™ ± Xq)) diverges log- 
arithmically at xq ~ > 7r/2; it should be cut off when 
Sx n ~ 1, i.e., at n/2 — xo ~ (LTsia((j>/2)/hv x ) . There- 
fore, within the logarithmic accuracy, 



(tanxo-F(7rn ± xq)) = 



In 
f 



hi)' 



2D 



LT 



(2n±l) 2 



tanh 



nhv x (2n ± 1) 
8LT 



For n = the integral diverges logarithmically at xo — > 
7r/2. At the lower limit, xq — > 0, the function F(x)tanx 
under the integral should be replaced with that contain- 
ing the resonance level taken from Eq. (|43l) . Therefore, 
we find 



(tan x F(x )) 



- ] b, 

7T 



LT J 8LT 7T w 



where 
B{0) = 



1 hv x x 4 nhv* 
tann =r tann ■ 



4LT ir 2 



8LT 

■° dx 



tan xq dxo 
Hv x x 



tanh 

o x 4LT 



Here x = yj x\ + (2LT/hv x ) 2 cos 2 ((/>/2) and LT/Hv 2 d < 
ip < 1. For low temperatures, 



In 



ho 



2D 



In(^), r|cos(^/2)|«T« 



T<r|cos(</)/2)| 




FIG. 4: Spectrum of bound states in a Josephson junction as a 
function of superconducting phase difference. Here we put (a) 
T/A = O.f, 2LA/hv F = 10, hk F v F /A = 10, DA/hv F = 1, 
= 0; (b) T/A = 0.1, 2LA/hv F = 30, hk F v F /A = 10, 
DA/hv F = 1,0 = 



It is independent of within the logarithmic accuracy. 
Performing the summation we find for T <C v 2 d / L 



i- 11 ( Lr2 1 r 

eR \hv 2 D 



(46) 



The current-phase relation is sinusoidal; this corresponds 
to the limit of low transparency junctions Z> 1 realized 
in short- lead contacts with LT/Hv 2 r> <C 1. The supercur- 
rent increases with lowering the temperature and satu- 
rates for T < r. For temperatures higher than Hv 2 d/L, 
the logarithm decreases to a value of order unity, and our 
approximation breaks down. Note that the current-phase 
relation, Eq. (|46l) , is similar to that obtained for a double 
barrier SINIS structure with resonant transmission (see, 
e.g., Ref.HJ. 

If the leads are very short, L <C fe^u/A, there exists 
only one level, Eq. (|4"2"j) . The estimates show that, within 
the logarithmic approximation, the current is dominated 
by the spectrum close to the resonance, Eq. (|4"5|) . and 
still has the form of Eq. (gBJ). 



2. Long leads 

Here we again restrict ourselves to the limit T <C A. 
For the contact with long leads, LT/hv x ^> 1 the spec- 
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FIG. 5: Junction made of the superconducting electrode 
placed on top of the 2D electron gas. 



trum consists of one level with eo < T satisfying Eqs. 
(|35|) . (|37|) and a series of levels with e n > T satisfying 
Eqs. (|35|). (p6| . Typical spectrum of bound states for a 
contact with long leads is plotted in Fig. @|b). 
For e < T we find from Eqs. (J35J) , ([57)1 



(2LT/hv x ) sin(0/2) sinhx 
V sin 2 a' + sinh 2 a; 



(47) 



in Fig. [5j The junction consists of a semi-infinite normal 
2D layer a part of which (with a length L) is covered by 
a bulk superconducting lead. The current through such 
contact can be written in terms of reflection coefficients 
obtained in Sec. IIIIBI Using the results of Ref. 17 
have 



</<: / % — [1-M 2 + M 2 ] 

v x >0 2 V 2 D 



1 

x[f (e-eV)-f (e)] 



where /o(e) = [e e / T + 1] 1 is the Fermi function. For 
T <C r the differential conductance is 

dJ 1 W 2 

— = — jf cos9[l-\r N \ 2 + \r A \ 2 } d6, (50) 



where x = (2LT/hv x )^l - e 2 /T 2 . If LT/hv x > 1 we where e = eV. 
have sinha; ^ 1 thus x — (2LT/hv x )s'm((f>/2), i.e., this 
energy state is as in a ballistic junction with a gap T, 

e 2 = r 2 cos 2 (0/2) . 

This level gives the usual expression for the supercurrent 
as in a ballistic contact with an induced gap T— , 



7rrsin(0/2) , rcos(0/2) 
I = tanh ■ 



ei? 



2T 



(48) 



One also needs to consider the contribution of levels 
with energies above the induced gap. For LT/hv x ^> 1 
Eqs. ([55)1 and yield x' = irn + 5x n where x' = 
(2LT/hv x ) v / e 2 /T 2 ~ I and 

nn + 8x n = b\ sm(5x n )\ . (49) 

Here we denote b 2 = (2LT/hv x ) 2 sin 2 (0/2)/ sin 2 a'. Ex- 
panding in b^ 1 <C 1 we find Sx n — ±mi/b — nn/b 2 . For 
e>f we have 



Eden 

n 



hv x \ - 
"2L ^ 



d 1 



The terms in Sx proportional to b^ 1 disappear. The num- 
ber of terms in the sum is N ~ b, being determined by 
the condition Sx < 1. Therefore, the sum over the states 
with e > r is of the order of (hv x /2L) <C T, and can thus 
be neglected. As a result, Eq. (|4*5|) is the full expres- 
sion for supercurrent through a contact with long leads. 
In this sense, the contact with long leads behaves as a 
contact without normal reflection, |rjy| = 0, when the 
bound states would only exist for e < e g . This limit was 
considered in Ref J^. 



A. Long leads. Andreev reflection 

Consider the case L — > oo such that the length is longer 
than the electronic mean free path, L»f. In this limit 
one can neglect the effect of the dead end and assume 
that Z = 0. 

For Z — we have zero normal reflection amplitude 



and the Andreev reflection coefficient 



r A 



2 /\u\ 2 which is independent of p y . For e > A the 



Andreev reflection coefficients becomes 



<( 1 + 7^)-f 2 ( 1 + 7&) 



r 2 



el 



For e < A we have 



e 2 1 



2/r 



-T 2 



el 



-We 2 1 



2F 



-V 2 



( 1 + -7A^) + V £2 ( 1+ vfc) 



r 2 



If e 2 < e 2 , one has \u\ = \v\ and \r A \ 2 = 1. If e -> A — 
we also have Ir^l 2 -> I . If T < A then for e g < e < A 
we have \r A \ 2 = r 2 /4e 2 . Since the Andreev reflection 
does not depend on 9, Eq. (0 yields for T < T 



dl 

Ro dV = 1 



r A 



IV. THE I-V CURVE IN A 2D/S JUNCTION 

The results obtained in the previous section can be 
used to describe the transport in a 2D/S junction shown 



where e = eV. The Andreev reflection coefficient for a 
long contact is shown in Fig. [75] The limit of a long con- 
tact with zero normal reflection was considered in Ref^i 
using a direct-contact model. 
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B. Finite-length leads. Oscillations of the 
resistance 

Here we discuss contacts with finite-length leads in the 
weak coupling limit, T <C |A|. For e < |A| when quasi- 
particles cannot escape into the leads, Eq. (|29[) gives 
1 - |rjv| 2 + \ta\ 2 = 2\r A \ 2 - Equation §§0§ yields 



Ro 



dl 

dV 



tt/2 



cos 6\r a {e)\ 2 dO . 



In the low- voltage region, e < e g , we have u* = v while 
Z is real. Therefore, 



FA 



4wV sinh 2 [2LVr 2 -e 2 /hv x ] 



Au 2 v 2 cosh 2 [2lVF 2 - e 2 /hv x ] - 1 ' 
In the weak coupling case e g = V and 

r}! = eT/A , 7/2 = T , Au 2 v 2 = T 2 /e 2 , e = e 
The differential conductance becomes 



R dI 



2 x 



4L 2 hv 2D 

75— In ' L ^/ hv 2D < 1 

I , LT/hv 2D > I 



It does not depend on voltage. For short leads LT/v2D <C 
1, the conductance is determined by the incident angles 
9 close to 7r/2 where the Andreev reflection is of the or- 
der of unity. For long leads LT/hv2D 3> 1 the Andreev 
reflection is complete as in the previous Section IIV Al 
For energies e g <C e < |A|, the Andreev coefficient is 



2uV[l-cos[2($ + -5-)}} 
+ v*- 2u 2 v 2 cos[2(<5+ - Jjj] 



If the energy is not very close to the superconducting 
gap, \e — A| » T 2 /A, we have v 2 -C u 2 , while u 2 = I. 
Therefore 



2v 2 [l-cos[2(<5 + -5_)] 



0.5 



"•' ; P-«.5\0.K , 



s/A 



0.5 



1.5 



and v 2 (e) = T 2 /Ae 2 . Therefore, for e g < eV < A, 



i?0^r = ts-X 



21n-L. 
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^ 2 



f 



7T . / - 7T 

-^sin(2V 

2y v 4 



V < 1 

7> 1 



where = 2LeV/tvV2D- The conductance exhibits os- 
cillations as a function of the bias voltage due to the 
geometrical quantization. 

Consider now the energies eV > A 3> T. Not very 
close to the gap edge A the Andreev reflection is small, 
while 



|rjv| 2 = exp[— 21m (<5+ — 5-)] = cxp 



ALT 



hv x Ve 2 - A 2 
For short contacts LT/hv2D C 1 we find 



dl 2ttLT 



dV hv 2 D Ve 2 - A 2 ' 



(51) 



Region close to the gap edge, |e — A| — > needs a special 
consideration. Here w 2 — > w 2 , so that the normal reflec- 
tion vanishes while the Andreev reflection grows up to 
\ta\ 2 — 1, and the conductance has a sharp peak as in 
the case of long contacts considered in Sec. IIV Al The 
differential conductance is plotted in Fig. [7J 

As was mentioned in Sec. 1111 Bl short leads are equiva- 
lent to the limit of low transmission tunnel contact. The 
differential conductance Eq. (|5T1) at voltages eV > A is 
proportional to the superconducting density of states, as 
it is usually the case for tunnel contacts. However, in 
addition to the peak at the superconducting gap edge, 
the differential conductance has a low-energy peak which 
characterizes the induced superconducting gap e g in the 
2D layer. In the voltage interval e g < eV < A, the 
conductance oscillates as a function of voltage due to 
the geometrical quantization in the regions where the 
leads overlap with the 2D layer. The contacts with long 
leads, LT/hv2D 3> 1, for e > A have an exponentially 
small normal reflection due to almost complete escape of 
quasiparticles into the superconductors, and thus their 
conductance coincides with that obtained in Sec. IIV Al 




FIG. 6: Andreev reflection coefficient as a function of energy 
for zero normal reflection, Z — 0. Here x = e/|A|. 



FIG. 7: Differential conductance G = Ro(dI / dV) as a func- 
tion of V/IOOA. Here T/A = 0.1 and 2LA/hv 2D = 10. 
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V. DISCUSSION 

We formulate the approach which can be used for spa- 
tially inhomogeneous and/or time-dependent problems 
associated with the induced superconductivity in low di- 
mensional electronic systems including a 2D electron gas, 
graphene layer, etc. This approach is based on the so- 
called Fano- Anderson model which describes the decay 
of a resonance state coupled to a continuum 1 ^. We con- 
sider a 2D system placed in a contact with a bulk su- 
perconductor and use the tunnel approximation for cou- 
pling between the superconductor and the 2D electron 
layer in a way similar to that used in Refsi 12 ' 13 for im- 
purities in a superconductor. We consider two particular 
examples of junctions made of a ballistic 2D electron gas 
placed under the superconducting electrodes of a finite 
length. For the case of a short symmetric S/2D/S junc- 
tion we find the bound states localized in the junction 
and calculate the supercurrent as a function of the lead 
length. Next we consider the 2D/S junction and find 



the IV curve for various lead lengths. We show that the 
differential conductance as a function of the bias volt- 
age exhibits peaks which correspond to the induced and 
bulk superconducting gaps. The differential conductance 
also shows oscillations as a function of voltage due to ge- 
ometrical quantization in the regions with the induced 
superconductivity. 
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